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ABSTRACT 

We present analysis of the asteroid surface density distribution of main 
belt asteroids (mean perihelion A ~ 2.404 AU) in five ecliptic latitude 
fields, — 17>/3(°)< + 15, derived from deep Large Binocular Telescope (LBT) 
V— band (85% completeness limit V = 21.3 mag) and Spitzer Space Telescope 
IRAC 8.0 /im (80% completeness limit ~ 103 yuJy) fields enabling us to probe 
the 0.5-1.0 km diameter asteroid population. We discovered 58 new asteroids in 
the optical survey as well as 41 new bodies in the Spitzer fields. The derived 
power law slopes of the number of asteroids per square degree are similar within 
each ~ 5° ecliptic latitude bin with a mean value of —0.111 ± 0.077. For the 
23 known asteroids detected in all four IRAC channels mean albedos range from 
0.24 ± 0.07 to 0.10 ± 0.05. No low albedo asteroids (p v < 0.1) were detected 
in the Spitzer FLS fields, whereas in the SWIRE fields they are frequent. The 
SWIRE data clearly samples asteroids in the middle and outer belts providing 
the first estimates of these km-sized asteroids' albedos. Our observed asteroid 
number densities at optical wavelengths are generally consistent with those de- 
rived from the Standard Asteroid Model within the ecliptic plane. However, we 
find an over density at (3 ^ 5° in our optical fields, while the infrared number 
densities are under dense by factors of 2 to 3 at all ecliptic latitudes. 
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INTRODUCTION 



The present main belt asteroid size distribution provides important constraints on mod- 
els of the original size distribution of planetesimals and their collisional evolution. The size- 
frequency distribution, main belt asteroid scale height, and albedo variations extant in a 
large, well-sampled population are also essential to model the size distribution of the Near 
Earth Asteroid (NEA) population and to determine the evolutio n of main belt obje cts as 



compared to those bodies injected into Earth crossing orbits (e.g.. iTrilling et al.l 120081 ) . 



Most main belt asteroids are found between 2.2 and 3.4 AU from the sun and at ecliptic 
latitudes nominally less than 20°. The size- frequency distribution of main belt asteroids with 
diameters < 1 km is essentially unknown due to observational flux limits. For asteroids with 
diameters > 1 km, the derived size-frequency distributio n estimates range from 7 x 10 5 km 



sized objects (optical Sloan Digital Sky Survey, SDSS; lAbazajian et al.l 12004 ; llvezic et al. 
20011 ) to (1.2 ± 0. 5) x 10 6 km-sized objec ts (derived from surveys with t he Infrared Space 



Obser vatory, ISO; iTedesco fe Desertll2002l ) within the main belt. Recently, IStapelfeldt et al. 



( 120061 ) have derived a size-frequency estimate for the > 1 km-sized bodies for two fields at 
0° and 10° ecliptic latitude from 24 /im photometry obta ined with the NASA Spitzer Space 
Telescope (Spitzer. iGehrz et al.l 120071 ; IWerner et al.l 120041 ) that concurs with values derived 
from analysis of the Sloan fields. While the infrared (IR) and optical methodologies yield 
comparable size estimates, these surveys are potentially biased as they are unable to detect 
(either through intrinsic sensitivity limits or survey design) fainter objects within the main 
belt. Thus, the size- frequency distribution may not be complete due to the preferential 
selection of a cross section of the asteroid population with high albedo values or large sizes. 



Bottke et al.l (120051 ) have modeled the collisional evolution of the size distribution of the 



main belt using the current size-frequency distributions within the main belt as well as the 
size-frequency distribution of NEAs to constrain the primordial asteroid distribution. Their 
best-fit models, based on the current size-frequency dist ribution, predict ~ 10 6 asteroids 



within the 0.3 to 1 km size range. Collisional models by lO'Brien fe Greenbergl (120051 ) pre- 
dict ~ 10 5 bodies in a similar size range within the main belt. However, the validity of the 
population extrapolation is based on an accurate understanding of the current size-frequency 
distribution necessitating deep observations to provide accurate number counts at both op- 
tical wavelengths (higher albe do objects) and IR (larger dark asteroids) wavelengths. The 
collisional evolution model by iBottke et al.l (120051 ) suggests that the primordial main belt 
contained ~ 150 to 250 times the current population of objects with diameters, D < 1000 km. 
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The IBottke et al.l (120051 ) model produces a parent population of si milar magnitude to that 
required by interpretation of the meteoritic record. For instance IWetherilll (119891 ) argues 
that the surface density within the main belt must be at least 100 times higher than that 
currently observed to be consistent with the meteoritic record. The total mass of known 
asteroids in the main belt is 5 x 10~ 4 M©. However, both the meteoritic record and colli- 
sional models require a primordial main belt mass of 7.5 x 10~ 2 M e to 1.25 x 10 _1 M®. The 
significant discrepancy between the mass of the primordial main belt and the current mass, 
primarily derived from large (D > 10 km) sized asteroids, may be resolved by extending 
the asteroid number census to properly encompass smaller sizes while including populations 
found at higher ecliptic latitudes. 

Here we present new observational results to constrain the size-frequency distribution 
of main belt asteroids and their ecliptic scale-height derived from a deep optical survey of 
select ecliptic fields obtained with the Large Binocular Telescope (LBT) as well as near- 
and mid-IR fields observed with Spitzer drawn from the data archive. Section [2] describes 
our observations and reduction techniques, $3] discusses the newly discovered asteroids and 
our derivation albedo, size-frequency distribution, and scale-height, while §H summarizes our 
conclusions. 



2. OBSERVATIONS & ARCHIVAL ANALYSIS 



2.1. LBT Observations 



Optical observations were obtained at the Large Binocular Telescope (LBT; iHill et al. 



20061 ) facility of the Mt. Graha m International Observatory with the blue c hannel of the 
Large Binocular Camera (LBC; Ragazzoni et al. 2006 ; Giallongo et al. 20081 ) and a single 
8.4 m mirror on various nights during 2007 January 16 through 24 UT as part of Science 
Demonstration Time (SDT) activities. The LBC is a wide-field imager incorporating four 
2048 x 4608 pixel CCD detectors with a 23' x 23' field of view (FOV) and a 0.23" per pixel 
platescale. Four fields at ecliptic latitudes of 0°, 5°, 10°, and 15° were observed using a series 
of 4 min exposures in the V-band (A = 0.55 /im; AA = 0.094 /im) under non photometric 
conditions with seeing between ~ 1.3"— 3.0". Fields were selected to be at solar elongations 
near 160°, such that sufficient asteroid motion on the plane of the sky could be detected 
within the ~ 3 hr observational baseline. Observations of each field consisted of 2 pointings, 
with one pointing offset by ~ 7.85', resulting in a total areal coverage per field of 0.2835 
square degrees. Table [Q provides complete observational details. 
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The data were reduced using standard IRAF0 routines in the mscred package. Images 
were trimmed, bias-subtracted, and flat-fielded using median sky flats created with twilight 
flats and additional data obtained during 2007 January SDT time. Astrometric solutions 
for each data frame were generated using IRAF routines, in conjunction with the stellar 



positions obtained from the USNO A2.0 astrometric catalog (IZacharias et al.ll2004l ). The 
solutions had a positional accuracy of 0.25" ± 0.08" for each field and the relative position 
of detected asteroids are referenced to this grid. 

Asteroid detection was performed using a thre e color method. Imag es with astrometric 



solutions were displayed in image display tool DS9 (I Jove fc Mandelll2003l ). with each individ- 
ual epoch (Table [I]) loaded into a different color table. Asteroids were identified by pairs or 
triplets of individually colored sources. From these identifications, absolute pixel coordinates 
for each target were obtained, which were subsequently used for photometric measurements. 
Due to the conditions under which the data was obtained, a high signal to noise of 10 is 
required for any asteroid detection. Secondary confirmation of this asteroid detection tech- 
nique was performed by registering individual epoch frames by the world coordinate system 
(WCS) and then subtracting, leaving positive and negative asteroid pairs. While this latter 
technique would be preferred over the three color technique, the seeing varied by 1.3" to 3.0" 
over each night causing the WCS subtraction technique to be less reliable. 

Photometry was performed using a circular aperture with a radius of 10 pixels for each 
asteroid coordinate. The photometry was calibrated from stellar sources in the fields using 
published V— band magnitudes taken from the Tycho and USNO A2.0/YB6 catalogs (mean 
photometric errors of the order < 0.1 mag). Forty stars (10 stars per chip) were identified 
by their astrometric positions and were used to find a photometric offset between the mea- 
sured and reported absolute values from Tycho/USNO A2.0. This zero-point photometric 
offset, which includes the mean photometric error, was then used to calibrate the asteroid 
photometry in each pointing as described in Tabled], column [7]. 

We detected 62 asteroids in the LBC V— band images, of which only four are previously 
known objects with orbital determinations. The 58 newly discovered asteroids detected in 
this survey were observed in either two or three epochs. Due to the field overlap between 
pointings, some asteroids were detected six times, allowing for both precise astrometry and 
photometry. Table [2] summarizes the number of asteroids found in each field as well as the 
extrapolated number counts per square degree at each latitude and the limiting magnitudes 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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for each field. The extrapolated optical number counts were obtained by multiplying the 
asteroid number counts in a field by the number 0.2835 square degree fields necessary to tile 
a 1 square degree field assuming a uniform asteroid density at each latitude. 



2.2. Spitzer Archival Fields 



Asteroid photometry and number counts in select ecliptic fields discussed in §3] were 
derived from new ground-based optical observations a nd point-sources e xtractions obtained 
from fields observed with the Infrared Ca mera (IRAC; lFaziq et al.ll2004T) and the Multiband 
Infrared Photometer for Spitzer (MIPS; iRieke et al.l 12004 ) as part of various IR survey 
programs retrieved in the Spitzer public archive. 

Post-pipeline (pipeline vl5) basic calibrated data (BCDs) of selected Spitzer IRAC 
fields were downloaded from the public archive and utilized for asteroid detection. Fields 
were selected from the ecliptic plane component of t he Spitzer First Look Survey (FLS, 
prog ram identification (P IP) 98; iMeadows et all 120041 ) . and SWIRE XMM-LSS fields (PID 
181; lLonsdale et aPl2003h . IRAC 8.0 fim data is preferable for asteroid detection due to the 
detector platescale (~ 1.2") and detector sensitivity. The 5a detection limit for a 500 sec 
IRAC 8.0 jum observation is 27 /xJy, enabling detection of small main belt asteroids with high 
signal-to-noise ratios (SNR) in short integration times. For example, the 8.0 fim flux from a 
1 km diameter main belt asteroid radiating as a blackbody with an orbital semi-major axis 
of 2.5 AU (3.2 AU) viewed at opposition is « 1096 /iJy (~ 20.1 /iJy). Once an asteroid 
candidate is identified at 8.0 [im its sky coordinate can then be used to examine pixels at 
the same location in the IRAC 3.6, 4.5, and 5.8 fim images. Although the flux densities for 
small-sized asteroids near the outer belt edge are not easily detectable with IRAC in single 
frames at shorter wavelengths, their motion will separate them from the confusion caused 
by faint extragalactic sources or point spread function (PSF) smearing of the telescope. 

Due to limits on observation durations (driven by on-board data storage issues and 
downlink frequency), the Spitzer deep large area surveys must consist of multiple epochs of 
data with multiple pointings. For asteroid detections this is advantageous as the observation 
duration limit for all IRAC observations is 6 hrs which allows observations to be repeated 
with a cadence of a few hours. Assuming circular Keplerian orbits within the asteroid belt, 
the angular velocity u (rad s _1 ), or rate of motion on the sky of an object near opposition is 



(1) 
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where G = 6.67 x KT 11 m 3 kg" 1 s" 2 , M is the mass of the sun, and R is the orbital distance 
in meters. For inner main belt asteroids with a semi-major axis of 2.2 AU, the expected rate 
of motion is 45" per hr, while the expected rate of motion for outer main belt asteroids is 
24" per hr. However, the actual rates and direction of motion on the sky may significantly 
differ from these values (e.g., near quadrature, when the FLS fields were observed at solar 
elongation ~ 115°). Nevertheless, movement on the sky of many tens of arcsec per hr enables 
one to detect asteroids in the field on timescales of hours. 

The FLS consisted of IRAC observations of two 0.13 square degree fields (10' x 48') 
centered on a solar elongation of 115° as seen from Spitzer, at ecliptic latitudes (3 = 0° and 
+5° on 2004 January 21 UT. In J2000 coordinates, this corresponds to field centers of R.A. 
= 12 h 03 m 22.03 s , decl. = -00°21'53.8" for the (3 = 0° field, and RA = 12 h ll m 20.10 s , decl. 
= +04°13'18.7" for the j3 = +5° field as viewed by Spitzer. To detect asteroid motion, each 
IRAC field was observed at three epochs separated by 70 min. 

For each of the latitude fields, we colored each epoch and coadded the IRAC 8.0 jum 
mosaics of the three epochs using a fixed stellar WCS to create a composite red, green, blue 
(RGB) image, which was then examined. In this RGB image, fixed sources appear white, 
and moving targets were seen as red-green-blue quasi-linear sequences of sources. The typical 
SNR for the detected sources was > 10. Using the RGB technique with the IRAC 8.0 //m 
data, we were able to identify 18 sources in the (3 — 0° field, and 16 sources in the (3 = +5° 
field. The faintest of these sources have 8. Aim fluxes of ^ 80 /iJ y and were detected with a 



SNR ~ 10 in both fields. For comparison, iMeadows et al.l (120041 ) identified 15 known (six of 



which we also recover) and 19 unknown asteroids their analysis of the FLS fields, reporting 
8.0 iim fluxes ranging from 5.6 to 0.09 mJy 

The SWIRE program wa s designed to explore large scale structure in the universe with 



the IRAC and MIPS cameras (ILonsdale et al.ll2004l ). This survey consisted of 6 fields selected 



to overlap those observed at depth in by the Hubble Space Telescope and the Chandra X-Ray 
Observatory. Two epochs of observations were obtained for each of the six fields to preen 
moving targets from final mosaics enabling high fidelity point source catalogs to be generated. 
For the purposes of our asteroid survey we used IRAC data obtained between 2004 July 23 
and 2004 July 28 in the SWIRE XMM-LSS field, located at (3 = -17° and encompassing an 
area of 9.1 square degrees. Due to the large area covered by this Spitzer survey and duration 
limitations of single observations, the IRAC data consists of 16 individual pointings of 0.44 
square degrees in two epochs. The two epochs of each pointing are offset by 3.75 hrs, thus 
significant asteroid motion can be detected. 

Asteroid detection in the SWIRE data sets was also performed by coloring and coadding 
IRAC 8.0 Atm mosaics of each epoch to create a composite red-blue (RB) image. In this RB 
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image, fixed sources appear purple, and moving targets were seen as red-blue quasi-linear 
sequences of sources. We identify a total of 46 sources in the field. The faintest sources have 
8.0 //m fluxes of ~ 300 fxJy and were detected with a SNR ~ 20 in both fields. Of the 46 
sources detected, only 14 asteroids have ground based detections and orbital elements. 

Fluxes in all four IRAC channels were obtained for all asteroids in the FLS and SWIRE 
datasets. Pixel coordinates for each asteroid were recorded at the time of detection and were 
used for photometry. Though the IRAC cameras have an offset between the 3.6/5.8 //m and 
the 4.5/8.0 jum FOVs, pipeline processed mosaics in each channel are aligned such that a 
given sky coordinate corresponds to the same pixel coordinate in all four channels. This 
allows use of one set of pixel coordinates for all sources in all 4 channels. Source photometry 
was performed using the IRAF task phot and a three-pixel radius circular aperture to measure 
the asteroid flux and a concentric four-pixel wide sky annulus to determine the median sky 
contribution to the total flux measured in the integrated object aperture. Final fluxes and 
flux uncert ainties were derived using th e aperture corrections contained in the IRAC Data 



Handbook (ISpitzer Science Centerll2006l ). Due to the offset in the IRAC channel FOVs, some 
asteroids are only detected in the 4.6 and 8.0 jttm data. The number of asteroids detected in 
all four IRAC bands is 40 in the SWIRE data and 24 in the FLS fields. The mean (median) 
8.0 asteroid fluxes are 315 ±78 /iJy (199 ±78 /iJy), 1276 ±680 fiJj (302 ±608 fiJy), and 
2448 ± 53 /iJy (1355 ± 53 ^xJy) for Spitzer FLS 0°, FLS 5°, and SWIRE fields respectively. 
The total number of asteroids detected in the Spitzer fields are summarized in Table [2j 

Detailed physical parameters for the newly discovered asteroids detected at optical (58) 
and IR (41) wavelengths are presented in on-line machine-readable catalogs. Illustrative 
catalog layouts are shown in Table [3] and Table HI The optical catalog, Table El is similar 
to the standard Minor Planet Center submission formai@. In our Table [3] the corresponding 
column format is used: [1] the asteroid name, [2] a true/false=l/0 flag indicating if asteroid 
is a new discovery, [3]-[5] are respectively the year, month and UT date of observation, 
[6]-[ll] are the right ascension and declination (J2000.0) of the detected target, [12] and 
[13] are the observed asteroid magnitude and associated uncertainty, while [14] is the filter 
used for the observations. The IR machine-readable catalog, Table HI is similar in column 
format where: [1] is the asteroid name, [2] is the Spitzer Program ID number, [3] is the 
astronomical observation request (AOR) key, [4]- [9] is the year, UT date and time at the 
start of an observation, [5] and [6] are the right ascension and declination in degrees, [10]- 
[13] are the right ascension and declination and ecliptic latitude and longitude in sexigesimal 
degrees, while [14] through [21] are the observed fluxes and flux errors in /^Jy measured 



2 http : //www . cf a . harvard . edu/iau/inf o/OpticalObs . html 
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for each asteroid in the four IRAC channels. Note that the fluxes in the table are not 
color-corrected. 



3. DISCUSSION 

The distribution of asteroids within the main belt provides a means to ascertain the 
mass distribution of material within the main belt and provides constraints for theoretical 
models developed to describe planet formation from protoplanetary disks encompassing our 
own proto-sun, as well as those extant around other young stellar systems. By probing 
the asteroid distribution in reflected optical light and in the thermal IR in coordinated 
observational programs, a more complete picture of the number of aste roids and sizes and 



albedos can be obtained. Current models of solar system dynamics (e.g.. iBottke et al.ll2005l ) 
use the distribution of absolute magnitudes of asteroids combined with a single geometric 
albedo to estimate the size-frequency distribution of asteroids as observationally determined 
diameters and albedos of main belt asteroids only exist for ~ 2400 objects. However, new 
survey data now available from fields observed with the Spitzer can be used to determine 
asteroid albedos and diameters for sizes under 1 km, enabling a more robust determination 
of albedo trends and size as a function of latitude for a more representative population of 
objects (e.g., a broader sampling of the total distribution). Our study describes the first 
results for main belt asteroids with diameters of 0.5 km. 



3.1. Completeness Limits 

To test completeness of our optical and Spitzer data sets, we used the IRAF task 
mkobjects to create synthetic point sources in data frames. To assess completeness in the 
optical data, synthetic point sources were added with the same seeing values measured in 
each frame and with coordinates that varied to mimic asteroid motion. To model the Spitzer 
images, synthetic objects were randomly placed into an image with motions equivalent to the 
median rates of a body in Keplerian orbit at the corresponding ecliptic latitude. Synthetic 
asteroids were then recovered using the same RGB technique originally used to detected 
asteroids in the observational images. Our completeness values are not only sensitive to 
source fluxes, but also to motion. For valid completeness detection, an asteroid must be 
detected in at least 2 of the 3 epochs of data. 

In our optical data, the 85% moving target completeness limit in all fields is V = 21.3. 
Because the moving target completeness samples the completeness due to (1) a limiting 
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flux and (2) motion, a point source completeness test was also performed. This exercise 
shows that the point source completeness at V — 21.3 is 90%. For all asteroids counts at 
magnitudes greater than 21.3, the difference in the two completeness values of 5%, is applied 
by assuming that the number of asteroids not detected due to motion is the same at all 
magnitudes. These completeness values have been used to obtain corrected asteroid number 
counts which are reported in column [3] in Table [2j Figure [1] shows the number counts per 
square deg function of magnitude in each optical field. 

The Spitzer IRAC 8.0 /im analysis, under the constraints that the synthetic asteroid rate 
of motion was comparable to the median rates of motion in a field at a given ecliptic latitude 
combined with the requirement of a minimum of 2 detections per synthetic asteroid yields 
an 80% completeness limit for the FLS and SWIRE data for F 8 ^ m = 103 u Jy. Our derived 



compl eteness value differs somewhat from that reported for the FLS by Meadows et al. 



( 120041 ) who cite 90% completeness at 100 /^Jy based on fixed target detections. With the 
completeness limit at 103 /iJy, we calculate that the number of asteroids in the FLS 0° and 
5° fields are 148 ±12 and 133 ±12 per square degree respectively, and 5 ±2 per square degree 
for the —17° SWIRE field. Figure [2] are histogram plots of the surface density of asteroids 
as a function of 8.0 /iin flux in the Spitzer FLS and SWIRE fields. 



3.2. Bulk size-frequency distribution 



Comparison of these two major data sets in different wavelengths is complicated in part 
by the potential for a disparate probe of asteroid sizes. Using the completeness limits in the 
V^-band and at IRAC 8.0 /mi we can estimate the minimum radii of asteroids detectable in 
our survey from 



2hv 3 1 R Q Rp v nD 2 cos(a) 



reflected & _ x r 16ttA 2 



(2) 
and 



thermal 



2eD 2 hv 3 

"'Spitzer (' 



. 7T „ 7T 



-cos 2 {(f)) cos {9 — a) d6 < 



(3) 



g kT ast — \ 



where 



T, 



ast 



;i - a)s ( 



ThV^SB 



cos6) A (cos9)*. 



(4) 



-10- 



The terms F re f\ ecte d and F t hermai are the observed y-band and IRAC 8.0 //m band fluxes 
(W m -2 Hz) respectively. Variables in the equation for reflected flux are: v is the frequency 
of observation, in this case V-band, R Q and T are the radius (m) and temperature (K) 
of the Sun, r is the asteroid heliocentric distance in meters, D is the asteroid diameter in 
meters and A is the asteroid geocentric distance, p v is the geometric albedo of the asteroid 
in U-band, R is the relative reflectance of the asteroid as measured with respect to the 
U-band reflectance and a is the observed phase angle in radians. Other variables in the 
equation for thermal flux are: dspitzer is the asteroid- Spitzer separation (AU), e is the asteroid 
emissivity which is assumed to be 0.9, A is the Bond Albedo, S & is the solar constant, r% 
is the asteroid heliocentric distance in AU, and rj is the beaming parameter. Adopting 
mean values for the known asteroids in our datasets = 2.82 AU, 5 = 1.89 AU (optical) 
and Th = 2.68 AU, dspitzer — 2.28 AU, (Spitzer asteroids; Table EJ), assuming rj = 1 with 
p v = 0.15 and cos(a) = 0.99, the minimum asteroid radii detected in our data are ~ 1020 m 
in the optical and ~ 450 m in the IRAC 8.0 /mi channel. Although Spitzer is far more efficient 
than optical imaging surveys at detecting smaller asteroids, asteroids identified from these 
two data sets comprise a set of bodies with similar size ranges within the main belt. 

The IR and optical data can be used to obtain the power law size-frequency distribution 
at 3 similar latitudes. In both the optical and the mid-IR, flux can be used as a proxy for size 
in obtaining the size-frequency distribution assuming that the Bond albedo is well behaved 
as a function of wavelength. By adopting a power law distribution for the surface density 
distribution oc f~ a , we use the completeness corrected counts to measure the power law 
slopes summarized in Table UJ The derived optical and IR slopes are similar at all latitudes 
within the uncertainties, with a mean value of a = —0.111 ± 0.077 for the surface density 
distribution of asteroids. However, the formal errors suggest that the optical data provides a 
better constraint on a. Our derived a- value s (Table [7J for the Spitzer fields differs somewhat 



from those cited by iMeadows et al.l (120041) . The differenc es can be attributed to how the 



asteroid photometry was conducted; IMeadows et al.l (120041 ) employed PSF-fitting (using, at 



that time, a poorly determined IRAC PSF), whereas we resort to aperture photometry to 
measure the flux density. In addition, diff erences in comp l etenes s also lead to variances in our 



results as opposed to those presented in IMeadows et al.l (120041) . We cite 80% completeness 



at 100 /xJys for moving-objects, while IMeadows et al.l (120041 ) quote 90% completeness at the 
same flux level for point sources. 

As the derived slopes of the surface density distributions are similar at optical wave- 
lengths, we employed a Kolmogorov-Smirnov (KS) test using the mean rates of motion of 
the asteroids at 0° and 15° ecliptic latitude to discriminate whether these bodies detected 
in the optical fields are drawn from the same size and magnitude population distributions. 
Unlike a y 2 test, data is not binned in the KS test and the KS test can be used on data 
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sets of different lengths. Due to the low number of asteroids detected in the 10° and 15°, 
binning of the data was rejected due to the associated loss of information and arbitrariness 
of defining bin sizes. The derived probability is 96.88%, strongly suggesting that the two 
fields represent the same distance (rate) population. Similarly, a KS test to determine the 
likelihood that the asteroids represent the same size distribution yields only a 41.32% prob- 
ability that the objects at 15° represent the same size distribution as those detected at 0° 
ecliptic latitude. This statistical inference suggests that the size distribution of asteroids 
at 0° and 15° are dissimilar; however, the veracity of this conclusion is limited by small 
number detection statistics of our data. Table [8] summarizes our probability analysis of the 
magnitude distribution (e.g., asteroid diameters) as a function of ecliptic latitude compared 
to the 0° field. Lastly, although the KS test indicates that the 0° and 15° optical fields are 
likely not the same population, they do appear to follow a distribution in which large, thus 
bright, objects are rare, whereas small and thus faint, objects are more numerous. 



3.3. Comparison to models 



Our observed asteroid frequency distribution can be juxtaposed with current models 
describing the evolution of the asteroi d populations in the solar syst e m, su ch as the Statistical 
Asteroid Model (SAM) developed bv iTedesco. Celling, fc Zappalal koO$ \. The SAM uses a 
set of 8603 asteroids with absolute magnitudes less than 15.75 (asteroids with diameters of 

> 2.5 km) to determine a size- frequency distribution of asteroids which is assumed to be 
smooth to diameters of 1 km. This model also incorporates an albedo distribution derived 
from 15 dynamical families and 3 "background" populat ions based on 1980 asteroids w hich 
have d iameter and albedo det e rminations from M SX (ITedesco. Egan. fc Pried |2002| ) and 
IRAS (jVeeder fc Tedescol If 9921 ; ITedesco et al.l 120021 ). A number of asteroids are excluded 
from this model, including those with inclinations greater than 25°, those with eccentricities 

> 0.3. 

The asteroid number counts observed with both Spitzer and the LBT are less than those 
predicted by the SAM. For a limiting flux of 0.06 mJy at 8 /im the SAM predicts 430 ± 40 
asteroids at 0° latitude and 250±20 asteroids at 5° latitude. On e reason the SAM may overes- 
timate the number of asteroids is the assumed faint flux limit of lTedesco. Cellino. fc Zappala 
(120051 ) which corresponds to an effective diameter of 0.6 km. This flux limit is fainter than 
the completeness for both fields, thus the SAM over predicts the number of small asteroids 
which would be observed by Spitzer. The overestimate of the number of asteroids may also 
be due to a fundamental assumption of the SAM - the power law slope of the size-frequency 
distribution is continuous to diameters less than 1 km. This is not observed in the optical. 
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The size-frequency distribution appears to fall off around 3 km (IJedicke &: Metcalfe! Il998l ); 
however, the reasons for this fall off are not well understood. Two possible reasons for the 
fall off in observed data are: 1) the bright limiting magnitude of most asteroid surveys ex- 
cludes detections of asteroids smaller than 1 km in diameter, or 2) the asteroids smaller than 
1 km may not be solid, but may instead be easily disruptible piles of rubble. If the predicted 
number of asteroids is scaled from the SAM with a model that yields the number of asteroids 
with diameters > 0.6 km being 2.4 times the number of asteroids with diameters > 1.0 km, 
the SAM still predicts 179 ± 17 asteroids in the 0° latitude field and 104 ± 8 asteroids in 
the 5° latitude field. The estimate at the high latitude is nearly in agreement with our FLS 
results; however, the SAM still overestimates the number of asteroids in the 0° field. The 
latter discrepancy may due to asteroids with high albedos (> 0.25) which are not detected 
in the IRAC channels because they are too cool and thus their thermal fluxes are below the 
detection thresholds in this data. 

For the LBT data which has limiting diameters of ~ 2 km, the number counts are 
within the uncertainty of the SAM prediction. For the SDSS survey which has a point 
source completeness limit of V = 22.2, the SAM predicts 115 ± 10 asteroids at 0° latitude 
and 75 ± 5 asteroids at 5° latitude. This completeness limit is fainter than ours by 0.9 mag; 
however, the model values are nearly in agreement with the asteroid number counts at these 
two latitudes. Unfortunately, the SAM does not predict asteroid counts at higher latitudes 
{(3 > 5°) and the current version of the SAM only uses asteroids with inclinations less than 
20° to make estimates on number counts. A subset of our newly detected asteroids have large 
inclinations. In fact this limitation of the SAM complicates derivation of the size-frequency 
distribution and physical characteristics of well known and well characterize d targets such 



as Pa llas and asteroid dynamical groups such as Hungarias and Phocaeas (jCarvano et al. 
200lh . 



3.4. Asteroid Diameters, Albedos, and Colors 



Direct comparison of mean and median IRAC 8.0 /mi fluxes ( §2.21) of asteroids in the 
two FLS fields and the SWIRE field, reveals that the SWIRE asteroids are, on average 
across the population, brighter than those asteroids at 0° and 5° latitude detected in the 
FLS fields, by a factor of two. This variance in brightness could be explained as a difference 
in asteroid diameter assuming that the populations have comparable albedos. We explored 
this possibility by estimating the diameters and deriving albedos for all asteroids which were 
detected in all four IRAC bands u sing the Near Earth Asteroid Thermal Model (NEATM; 
Delbo fc Harris! 120021 ; iDelbol 120041 ) . In this model, the thermal flux is dependent upon the 
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sub-solar temperature and the temperature distribution of the surface of the asteroid. The 
latter is solely depende nt on the albedo a nd 77 as mo del var i ables. Thus, our implementation 
of NEATM follows the belbo & Harry (|2002h and belbol ((20oJ prescription whereby the 
model fits the IR fluxes and the optical absolute magnitude of asteroids (if) by varying the 
geometric albedo (through a ^-minimization) until a best fit is found. The Spitzer fluxes 
were first color-corrected acco rding to the prescription described in the IRAC Data Handbook 
( jSpitzer Science Centerl 120061 ) . resulting in the measured fluxes being divided by 1.1717 and 
1.1215 at 5.8 /xm and 8.0 fim , respe ctively. The asteroid diameter is then calculated by the 
relation of iFowler fc Chillemil (119921 ) which uses only the geometric albedo and the absolute 
magnitude as input values. Because IRAC asteroid spectral energy distributions are compos- 



ites o f both thermal emission and reflected solar light (e.g.. iMueller. Harris, fc Fitzsimmons 
20071 ). we derived albedo and diameter estimates of the asteroids using only data from IRAC 
channels 3 and 4. From the channel 3 and 4 fits, we find that the percentage of reflected solar 
light in our channel 1 photometry to be ~ 56% and ~ 14% in the channel 2 photometry. Our 
modeling results for known asteroids detected in our IRAC fields are summarized in Table [6j 
The mean diameters of asteroids in SWIRE is 5.91 km while the mean asteroid diameter in 
the FLS is 2.28 km. 

Analysis of the mean geometric albedos and derived sizes for these objects suggests 
they divide between the two populations. The mean geometric albedo for asteroids in the 
FLS fields is 0.24 ± 0.07 as opposed to 0.10 ± 0.05 for asteroids in the SWIRE fields. No 
low albedo asteroids (py < 0.1) were detected in the Spitzer FLS fields, whereas in the 
SWIRE fields they are frequent. Evidently the surface composition of the asteroids in the 
SWIRE fields contains more carbonaceous material than the asteroids in the FLS, as car 



bonaceous chondrites ha ve significantly lower albedos in the 0.03 to 0.11 range (lGaffeylll976 
Johnson fc Fanaldll973l ). 



However, the marked dichotomy in derived albedos is most likely a result of the two 
Spitzer surveys sampling different asteroid taxonomy classes within the general asteroid 
belt. An extensive analysis of ~ 8 8, 0000 objects in t he Sloan Digital Sky Survey Moving 
Object Catalog (SDSS MOC) 4 by iParker et al.l (120081 ) demonstrates that there is a strong 
correlation between SDSS color, asteroid taxonomy (C-, S-, V-type), and orbital elements. 
In Fig. [3] we plot the proper orbital elements a (orbital semi-major axis) vs. i (orbital 
inclination) of our SWIRE (green squares) and FLS (red triangles) detected asteroids overlaid 
on 207,942 n umbered asteroids (black dots) with orbital elements in the ASTORB files 
(lBowellll200ll ). The location of the major Kirkwood gaps as defined by IParker et al.l (120081 ) 
are indicated, and the colored symbols used to identify the Spitzer-detected asteroids are 
coded by derived albedo. Open symbols are for Spitzer asteroids with derived albedos, 
Pv^O.l, while filled symbols identify asteroids with py > 0.1. The majority (~ 80%) of 



-14- 



the low albedo asteroids detected in the Spitzer data reside at high inclination in the outer 
belt. Primitive asteroids, with C-type colors and generally low albedos (py < 0.1), are 
seen to dominate the asteroid population in the outer belt, while S-type asteroids with p v 
ranging from ~ 0.15 to 0.2 tend to be more dominant in the middle and inner asteroid belt 



( jParker et al.l 120081 ) . Our SWIRE data clearly samples asteroids in the middle and outer 
belts, and in fact, for the first time, we are able to estimate the albedos of km-sized asteroids 
in the outer belt (Table [6]). 



The low albedos (pv ^0.1) of the 
lar to the albedos deri ved for Jupiter-family comets, 0.02 



6 km-sized outer belt SWIRE asteroids are simi- 

< 



< 



0.06 ([Fernandez et al. 



20081 ; lLamy et al.ll2004l ). Recent modeling by lLevison et al.l (120081 ) investigating the capture 
likelihood of cometary planetesimals into the asteroid belt (leading the establishment of a 
D-type population) suggests that the inner edge of the D-type population is near a ~ 2.6 AU. 
Whether the population of low albedo asteroids discovered in our SWIRE fields discussed 
here are organic-rich, primitive objects can only be verified with follow-up spectroscopy. 

From flux measurements in all 4 IRAC channels, we can also begin to constrain the 
colors of asteroids. This is beneficial for both the understanding of bulk asteroid albedo 
variations as a function of size and other orbital parameters such as orbital semi-major axis 
or inclination and for the future exclusion of asteroids from galactic or extragalactic catalogs. 
Figure H] is an adaptation of the work by lLacy et al.l (120041 ) which shows the IRAC colors of 
fixed astronomical objects where we have now included asteroid color data derived from the 
our FLS and SWIRE field photometry. Asteroids mainly reside in one isolated locus in color 
space. The size of this locus is dependent upon both temperature and albedo. The shortest 
wavelength observations with IRAC at 3.6 iim trace the fraction of reflected incident solar 
flux, thus traces albedo, while the 8.0 /im to 4.5 /im ratio traces the thermal emission and 
thus the albedo and distance of an asteroid. 



4. CONCLUSION 

Our combined optical survey of ~ 0.96 square degrees and Spitzer IRAC 8.0 //m fields 
encompassing from 0.26 to ~ 9.1 square degrees of sky at five latitudes perpendicular to 
the ecliptic plane resulted in detection of 118 main belt asteroids, of which 91 are newly 
identified objects either at optical or infrared wavelengths. The optical and mid-IR asteroid 
counts and fluxes demonstrate that the slope of the size-frequency distribution is consistent 
with that of a similar size population at 0° and 15° ecliptic latitude. The derived power law 
slopes of the asteroid surface density distribution are similar at within each ~ 5° ecliptic 
latitude bin with a mean slope of —0.111 ±0.077. The observed asteroid number densities at 
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optical wavelengths are generally consistent with those derived from the Standard Asteroid 
Model within the ecliptic plane. However, we find an over density at (3 > 5° in our optical 
fields, while the infrared number densities are under dense by factors of 2 to 3 at all ecliptic 
latitudes. 

For the 28 known asteroids detected in all four IRAC channels, mean albedos range from 
0.24 ± 0.07 to 0.10 ± 0.05, statistically suggesting that these are two different populations. 
Our SWIRE data clearly samples low albedo asteroids in the middle and outer belts, and in 
fact, for the first time, we are able to estimate the albedos of km-sized asteroids in the outer 
belt. 

In order to increase our number statistics and lower the errors, we are completing the 
data-mining of additional, extant multi-epoch survey Spitzer IRAC and MIPS dataset and 
are undertaking complimentary ground based optical observations with the LBT and the 
2.3-m Bok Telescope at Kitt Peak. These data will provide a survey area more than 3 
times larger than that described in this manuscript and will also expand analysis of asteroid 
number counts to ecliptic latitudes of 20°. 

The errors in the slope of the size-frequency distribution cannot be decreased by greater 
number statistics alone. The use of fluxes as proxies for asteroid sizes does not account for 
variations in albedo; therefore, additional mid-IR data is required to obtain albedos and 
sizes for a large number of asteroids. We in the process of completing the data-mining of 
additional, extant multi-epoch survey Spitzer IRAC and MIPS dataset while undertaking 
complimentary ground based optical observations with the LBT and the 2.3-m Bok Telescope 
at Kitt Peak in order to increase our number statistics and lower the completeness statistics 
uncertainties. These data will provide a ecliptic survey area more than 3 times larger than 
that described in this manuscript and will also expand analysis of asteroid number counts to 
ecliptic latitudes of 20°. These data will extend the number of asteroids with accurate sizes 
to diameters much smaller than the 10 km IRAS detection limit enabling reliable number 
counts for asteroids with diameters < 1 km and assessment of whether the fall off in the size- 
frequency distribution of main belt asteroids distribution with diameters < 1 km is a real 
signature, or an artifact introduced by the limiting magnitudes of current asteroids surveys. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under a con- 
tract with NASA. Support for this work was provided by NASA through an award issued by 
JPL/Caltech. Support for this work also was provided by NASA through contracts 1263741, 
1256406, and 1215746 issued by JPL/Caltech to the University of Minnesota. E.R. and 
C.E.W. also acknowledge support from the National Science Foundation grant AST-0706980. 



-16- 



This work is also based on data acquired using the Large Binocular Telescope (LBT). The 
LBT is an international collaboration among institutions in the United States, Italy and Ger- 
many. LBT Corporation partners are: The University of Arizona on behalf of the Arizona 
University system; Istituto Nazionale di Astrophysica, Italy; LBT Beteiligungsgesellschaft, 
Germany, representing the Max-Planck Society, the Astrophysical Institute Potsdam, and 
Heidelberg University; The Ohio State University, and The Research Corporation, on behalf 
of The University of Notre Dame, University of Minnesota, and University of Virginia. The 
authors also acknowledge the detailed and very helpful comments of anonymous referees 
whose insights greatly improved the manuscript. 

Facilities: Spitzer (IRAC, MIPS) LBTO (LBC 4048 CCD) 
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Table I. LBT/LBC OBSERVATIONAL SUMMARY 





Field 


Field 






Zero-point 


Pointing Name 


Center 


Center 


Date 


Time 


Seeing Offset 




RA (J2000) 


Dec (J2000) 


(2007 UT) 


(UT) 


(") (mag) 



degree field, pointing 1 



Epoch 1 


6:14:21 


68 


+23:25:15 


Epoch 2 


6:14:21 


68 


+23:25:15 


Epoch 3 


6:14:21 


68 


+23:25:15 


degree field, pointing 2 








Epoch 1 


6:14:22 


33 


+23:40:18 


Epoch 2 


6:14:22 


33 


+23:40:18 


Epoch 3 


6:14:22 


33 


+23:40:18 


5 degree field, pointing 1 








Epoch 1 


6:14:56 


16 


+28:25:17 


Epoch 2 


6:14:56 


16 


+28:25:17 


Epoch 3 


6:14:56 


16 


+28:25:17 


5 degree field, pointing 2 








Epoch 1 


6:14:56 


78 


+28:40:19 


Epoch 2 


6:14:56 


78 


+28:40:19 


Epoch 3 


6:14:56 


78 


+28:40:19 


10 degree field, pointingl 








Epoch 1 


6:15:33 


49 


+33:25:15 


Epoch 2 


6:15:33 


49 


+33:25:15 



49 


Jan 


16 


07:55:20 


2 


13 


7 


37+0 


10 


49 


Jan 


16 


08:00:01 


2 


74 


7 


50+0 


11 


49 


Jan 


16 


09:51:30 


2 


52 


7 


40+0 


13 


72 


Jan 


16 


08:01:21 


2 


19 


7 


68+0 


11 


72 


Jan 


16 


09:13:54 


2 


67 


7 


81+0 


06 


72 


Jan 


16 


09:57:22 


2 


67 


7 


78+0 


06 


94 


Jan 


25 


02:49:41 


1 


50 


7 


43+0 


04 


94 


Jan 


25 


03:49:01 


1 


64 


7 


47+0 


04 


94 


Jan 


25 


06:39:50 


2 


20 


7 


66+0 


04 


50 


Jan 


25 


02:55:22 


1 


50 


7 


48+0 


04 


50 


Jan 


25 


03:54:47 


1 


75 


7 


56+0 


04 


50 


Jan 


25 


06:45:37 


1 


71 


7 


66+0 


04 


14 


Jan 


24 


02:17:27 


1 


44 


7 


36+0 


03 


14 


Jan 


24 


04:34:49 


2 


12 


7 


62+0 


03 



Table 1 — Continued 





Field 


Field 








Zero-point 


Pointing Name 


Center 


Center 


Date 


Time 


Seeing 


Offset 




RA (J2000) 


Dec (J2000) 


(2007 UT) 


(UT) 


(") 


(mag) 


Epoch 3 


6:15:33.49 


+33:25:15.14 


Jan 24 


08:10:00 


2.64 


7.67+0.03 


10 degree field, pointing 2 














Epoch 1 


6:15:34.08 


+ 33:40:19.96 


Jan 24 


02:23:01 


1.44 


7.33+0.03 


Epoch 2 


6:15:34.08 


+ 33:40:19.96 


Jan 24 


04:40:57 


2.16 


7.55+0.03 


Epoch 3 


6:15:34.08 


+ 33:40:19.96 


Jan 24 


08:15:42 


2.65 


7.55+0.03 


15 degree field, pointing 1 














Epoch 1 


6:16:17.02 


+38:25:09.47 


Jan 24 


02:05:07 


1.35 


7.35+0.03 


Epoch 2 


6:16:17.02 


+38:25:09.47 


Jan 24 


04:13:06 


1.79 


7.60+0.02 


Epoch 3 


6:16:17.02 


+38:25:09.47 


Jan 24 


07:47:27 


2.62 


7.68+0.03 


15 degree field, pointing 2 














Epoch 1 


6:16:15.22 


+38:40:11.62 


Jan 24 


02:10:49 


1.28 


7.24+0.03 


Epoch 2 


6:16:15.22 


+38:40:11.62 


Jan 24 


04:19:06 


2.25 


7.60+0.03 


Epoch 3 


6:16:15.22 


+38:40:11.62 


Jan 24 


07:54:28 


3.08 


7.64+0.03 
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Table 2. SUMMARY OF ASTEROID DETECTIONS BY LATITUDE AND NUMBER 

COUNT STATITICS 



Field 


Number 


Number Total Number 


Limiting 


Ecliptic 


Asteroids 


Asteroids 


Asteroids 


Magnitude 


Latitude 


Detected 


Per Square Degree 8, Per Square Degree b 


or Flux 




LBT U-Band 


0° 


24 


85 


91± 10 


22.3 


5° 


21 


74 


80± 9 


22.0 


10° 


11 


39 


41± 6 


22.2 


15° 


6 


21 


23± 5 


22.2 


FLS IRAC 8/xm 


0° 


24 


85 


91± 10 


25.7 nJy 


5° 


21 


74 


80± 9 


25.7 /xJy 


SWIRE IRAC 8^m 


-17° 


11 


39 


41± 6 


35.4 fiJy 



a Directly calculated from the number of asteroids detected where each LBT 
field covers 0.24 square degrees of sky, the FLS fields are 0.13 square degrees, 
and the single SWIRE field is 9.1 square degrees in extent. 

b Number of asteroids including completeness. 



Table 3. OPTICALLY DETECTED ASTEROIDS - MACHINE READABLE TABLE a 



Asteroid 


Discovery 


Year 


Month 


Day 


RA a 


RA 


T\ A. 

RA 


Dec D 


Dec 


Dec 




Mag 




Name 


Flag 






(UT) 


(hrs) 


(min) 


(sec) 


(degr) 


(min) 


(sec) 


Mag 


Error 


Filter 


sdtOl 


1 


2007 


01 


16.33147 


06 


14 


03.68 


+23 


21 


03.0 


19.72 


0.03 


V 


sdtOl 


1 


2007 


01 


16.38195 


06 


14 


01.21 


+23 


21 


08.7 


20.34 


0.04 


V 


sdtOl 


1 


2007 


01 


16.41215 


06 


13 


59.70 


+23 


21 


12.1 


20.16 


0.04 


V 


sdt02 


1 


2007 


01 


16.33147 


06 


13 


56.23 


+23 


15 


27.2 


20.94 


0.13 


V 


sdt02 


1 


2007 


01 


16.38195 


06 


13 


53.82 


+23 


15 


26.8 


21.54 


0.10 


V 


sdt02 


1 


2007 


01 


16.41215 


06 


13 


52.35 


+23 


15 


27.0 


21.46 


0.04 


V 



a Full catalog available on-line. 

b RA and Declination given in J2000 coordinates. 



Table 4. SPITZER DETECTED ASTEROIDS - MACHINE READABLE TABLE a,b 



Asteroid AOR RA Dec Ecliptic Ecliptic 

Name PID Key YR MN Day Hi Min Sec J2000 J2000 Long Lat F 3 6 F| r ^" F 4 5 P| r 5 r F 5 8 F| r |" F g FfJ' 

(UT) (°) (°) (°) (°) (,xjy) (pjy) (pjy) ( M Jy) ( M Jy) ( M Jy) ( M Jy) (pjy) 



2001QY160 


98 


6095104 


2004 


01 


21 


02 


35 


25.728 


181.0084 


-0.0412 


180.9416 


0.3633 


89.4 


8.1 


18.9 


14.4 


95.7 


45.7 


939.9 


179.4 


Oh 


98 


6095104 


2004 


01 


21 


02 


35 


25.728 


181.0051 


-0.1183 


180.9693 


0.2913 


2.3 


2.2 


7.7 


6.9 


65.3 


35.5 


548.8 


132.8 


2004BS160 


98 


6095104 


2004 


01 


21 


02 


35 


25.728 


181.0004 


-0.1816 


180.9902 


0.2313 


3.7 


3.6 


1.7 


1.7 


52.5 


52.2 


103.2 


49.2 


2000AX136 


98 


6095104 


2004 


01 


21 


02 


35 


25.728 


180.9245 


-0.1815 


180.9204 


0.2012 


4.8 


4.6 


103.7 


47.2 


240.9 


79.9 


2526.7 


305.6 



a Full catalog available on-line. 

b Observed fluxes, no IRAC color-correction applied. 



I 

to 

I 



Table 5. OBSERVED ASTEROID IRAC FLUXES a 



Provisional 



Designation 


F 3 .e 


F4.5 


F5.8 


Fg.o 




(/*Jy) 


(A*Jy) 


(/*Jy) 


(A*Jy) 




SWIRE Asteroids 


2004TH222 


18. 5± 14.1 


53. 6± 30.6 


241. 0± 80.3 


1682. 6± 246.3 


2005YV181 


6.0± 5.6 


16.23± 12.8 


128.4± 55.2 


933.8± 179.0 


2005YV50 


16. 7± 14.2 


15. 7± 12.4 


64.5± 36.0 


570.6± 136.2 


1999WA1 


45. 1± 27.1 


113.4± 50.0 


938.5± 172.8 


7314.0± 532.0 


2004TO8 


109.2± 48.7 


259.0± 82.6 


1684.4± 237.5 


10301. 9± 634.3 


2006AP1 


5.8± 5.4 


32.8± 21.7 


263.8± 84.7 


1677.9± 245.7 


2005WX162 


10.4±8.9 


38.0± 24.1 


311. 5± 93.5 


2287.0± 290.30 


1999VY171 


75. 7± 38.5 


383.9± 103.8 


2867. 7± 314.2 


18500.5± 856.9 


1999JZ7 


4.0± 4.9 


25.5± 18.1 


196.8± 71.2 


1512.8± 232.3 


2000 VP 1 


15.9± 12.7 


25.0± 17.8 


113.9± 51.2 


775.1± 161.3 


2005YH117 


6.7± 6.2 


48.0± 28.4 


470.0± 118.0 


3570.0± 366.2 


2002FE15 


12.0± 10.1 


47.9± 28.4 


355.7± 100.8 


4308.9± 403.4 


2005YX74 


83. 1± 42.7 


118.5± 52.0 


538.9± 127.5 


4465.2± 411.6 


2002CN96 


15.5± 12.3 


65.0± 35.3 


584.1± 133.3 


3106.3± 343.6 



FLS Asteroids 



2001QY160 b 8.9± 7.8 31.8± 22.6 101. 3± 47.5 905.4± 175.8 



Table 5 — Continued 



Provisional 



Designation 


Hoc 

o.o 


I -1 A EC 


H r o 
- 1 - O.O 


Hon 
x o.U 




( z/Jv 1 


( 11 Jv J 


Kr^J ) 


( ziJv 1 


2000AX136 b 


114.8±50.4 


217.6± 74.2 


1228.2± 200.3 


8993. 1± 593.3 


2004BH160 


9.1± 9.3 


11.4± 10.9 


11. 5± 11.7 


73.6± 39.6 


2001QD49 b 


7.1± 7.4 


11. 2± 12.2 


16.1± 16.0 


109.9± 51.3 


2002WL7 b 


18.2± 13.9 


50.0± 29.1 


497.6± 121.4 


3239.4± 347.8 


1999FZ24 b 


41. 7± 25.9 


45.2± 27.3 


211. 0± 74.1 


2116.3± 277.6 


2001RP137 b 


8.3± 7.6 


14.3± 14.2 


34.9± 24.1 


273.5± 89.2 


2004EU91 


9.3± 9.2 


7.2± 6.8 


35.7± 26.5 


116.5± 53.3 


2004CA105 


61. 7± 35.8 


54.6± 33.7 


59.9± 34.8 


487.4± 124.5 


2002RG106 


12.2± 10.4 


5.3± 5.3 


22.9± 17.6 


228.3± 80.5 


2004BU99 


6.3± 5.8 


7.9± 7.4 


40. 3± 25.9 


398.9± 110.9 


2004BS160 


8.1± 8.6 


3.1± 3.1 


21.4± 30.7 


101. 9± 49.0 


2001SB182 


14.1± 11.9 


11. 8± 11.0 


32.6± 22.6 


184.0± 70.1 


2004FG8 


6.1± 7.2 


7.7± 7.0 


25.5± 19.2 


185.5± 70.7 



a Observed fluxes, no IRAC color-correction applied. 
b Asteroids also reported in iMeadows et all (120041 ) . 



Table 6. STANDARD THERMAL (NEATM) MODEL DERIVED ASTEROID DIAMETERS AND ALBEDOS 



Provisional 


Heliocentric 21 


Spitzer a 


Phase 


Absolute 13 


NEATM 


Geometric 


Beaming 


Model 


Designation 


Distance 


Distance 


Angle 


Magnitude 


Diameter 


Albedo 


parameter 


x 2 




(AU) 


(AU) 


(°) 


(H) 


(km) 










SWIRE Asteroids 


2004 1H222 


2.08 


1.79 


on a 


15.4 


2.83 


0.14 


1.69 


0.0002 


2005 YV181 


3.27 


3.09 


18.20 


15.5 


3.91 


0.07 


0.79 


0.0157 


2005 YV50 


2.38 


2.14 


25.36 


16.0 


2.45 


0.12 


1.66 


0.0087 


1999 WA1 


2.82 


2.60 


21.25 


12.7 


9.16 


0.18 


1.01 


0.0225 


2004 T08 


2.66 


2.44 


22.57 


14.5 


9.02 


0.03 


1.04 


0.5826 


2006 API 


2.91 


2.71 


20.53 


15.5 


3.59 


0.09 


0.77 


0.0035 


2005 WX162 


2.83 


2.61 


21.20 


15.3 


5.10 


0.05 


1.03 


0.0141 


1999 VY171 


2.56 


2.34 


23.48 


14.0 


10.92 


0.04 


1.04 


0.6644 


1999 JZ7 


2.33 


2.08 


26.02 


14.0 


3.39 


0.39 


1.35 


0.0308 


2000 VP1 


2.43 


2.21 


24.83 


15.7 


1.99 


0.23 


1.03 


0.0143 


2005 YH117 


2.91 


2.73 


20.57 


15.4 


6.73 


0.03 


1.01 


0.0161 


2002 FE15 


3.39 


3.25 


17.54 


14.5 


12.41 


0.02 


1.05 


0.3090 


2005 YX74 


2.85 


2.64 


20.99 


14.8 


5.35 


0.08 


1.05 


0.0859 


2002 CN96 


2.95 


2.74 


20.29 


14.8 


5.88 


0.06 


0.87 


0.9232 


FLS Asteroids 


2001 QY160 


2.96 


2.39 


17.74 


15.0 


3.26 


0.17 


1.02 


0.0241 



Table 6 — Continued 



Provisional 


Heliocentric 11 


bpitzer 


Phase 


A 1 l i K 

Absolute" 


NEATM Geometric 


Beaming 


TV T 11 

Model 


Designation 


Distance 


Distance 


A 1 

Angle 


Magnitude 


Diameter 


Albedo 


parameter 


2 

X 




(AU) 


/AT T\ 

(AU) 


( ) 


/TT\ 

(H) 


(km) 








2f)f)f) AX136 


2.40 


1 80 

X . VJ VJ 


22.02 


13 1 


6.89 


0.21 


1 33 

X . uu 


0002 

VJ . VJ VJ VJ Zj 


2004 BH160 

Zj VJ Ui iJlllUu 


2.42 


1.82 


21 86 

Zj X . O VJ 


18 4 


0.51 


0.30 


1 03 


0001 

VJ . VJ VJ VJ X 


2001 OD49 


2 77 

— ■ i I 


2.19 


18 98 


17 5 


0.94 


0.20 


1 00 

X . VJVJ 


0048 

VJ . VJ VJtiO 


2009 WL7 

iWuZi VV XJ ( 


? 69 

Zj . VJ C/ 


2 1 1 

Zj. 1 1 


19 55 


1 5 3 


4.48 


0.07 


1 05 


1406 


iqqq p7;24 


3 10 


2 53 


16 93 


13 8 


5.81 


0.16 


1 03 

X . V J (J 


1115 

VJ . X X X <LJ 


2001 RP137 


2.93 


2.36 


17.90 


16.7 


1.72 


0.13 


1.02 


0.0012 


2004 EU91 


2.40 


1.80 


22.04 


18.1 


0.51 


n 


0.76 


0.1778 


2004 CA105 


2.80 


2.22 


18.77 


14.7 


2.27 


0.44 


1.03 


0.0062 


2002 RG106 


2.69 


2.10 


19.58 


16.1 


1.30 


0.38 


1.045 


0.108 


9004 RTTQQ 


2.17 


1.54 


24.53 


17.4 


1.84 


0.06 


2.50 


0008 


2004BS160 


2.41 


1.81 


21.99 


18.6 


0.57 


0.20 


1.01 


0.0178 


2001SB182 


3.02 


2.46 


17.36 


15.6 


1.39 


0.53 


0.75 


0.0636 


2004FG8 


2.01 


1.37 


26.61 


18.5 


0.74 


0.13 


1.90 


0.0001 


a Heliocentric 


r h , and «5piteer-asteroid distances at epoch of observation. 








b Absolute 


Magnitudes 


are 


values 


given 


in the 


JPL Horizons 


database, 



http : //ssd . jpl . nasa . gov/?hor izons 
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Table 7. POWER LAW SLOPES OF THE SURFACE DENSITY DISTRIBUTIONS 







Ecliptic 


Slope 


Dataset 


Bandpass 


Latitude 


a 


LBT 


V 


0° 


-0.118 ±0.070 


LBT 


V 


5° 


-0.161 ±0.074 


LBT 


V 


15° 


-0.131 ±0.093 


FLS 


IRAC Bfjcm 


0° 


-0.083 ±0.145 


FLS 


IRAC Bfjcm 


5° 


-0.100 ±0.240 


SWIRE 


IRAC 


-17° 


-0.075 ± 0.340 
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Table 8. ASTEROID POPULATION PROBABILITIES, MEAN AND MEDIAN 

MAGNITUDES 



Field Ecliptic 



Latitude Probability 3, Median Mean 



(°) 


(%) 


(mag) 


(mag) 





100 


20.6 


20.45 


5 


66.24 


20.5 


20.41 


10 


38.89 


19.8 


20.00 


15 


41.32 


20.8 


20.65 
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LBT distribution ot 




20 

V Bond Magnitude 



LBT distribution ot 5 deg 



17 




20 

V Bond Mognitude 



17 



Fig. I. — Magnitude histogram distribution for LBT asteroids at 0° ecliptic latitude ([a], 
upper) and 5° ecliptic latitude ([b], lower). 
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LBT distribution at 10 deg 




21 20 

V Bond Magnitude 



LBT distribution ot 15 deg 



17 




21 20 

V Bond Mognitude 



17 



Fig. 1. — continued: Magnitude histogram distribution for LBT asteroids at 10° ecliptic 
latitude ([c], upper) and 15° ecliptic latitude ([d], lower). 
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FLS deg field 




-i 
Log Flux (mjy) 

FLS 5 deg field 




-1 

Log Flux (mjy) 

SWIRE -17 deg field 




-1 

Log Flux (mjy) 



Fig. 2. — Magnitude histogram distribution for Spitzer FLS asteroids at 0° ecliptic latitude 
(upper), at 5° ecliptic latitude (middle), and at 17° ecliptic latitude for the SWIRE field 
(bottom). 
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Fig. 3. — Proper orbital elements a (orbital semi-major axis) vs. i (orbital inclination) of 
asteroids. The Spitzer SWIRE and FLS asteroids are denoted respectively by green squares 
and red triangles, while 207,942 numbered asteroids with orbital elements in the ASTORB 
files are plotted with black dots. Open colored symbols indicate asteroids with derived 
albedos, pv less than or equal to 0.1, while filled sym bols denote a s teroid s with py > 0.1. 
The location of the major Kirkwood gaps as defined by lParker et al.l (120081 1 are indicated by 
the labels and vertical arrows. 
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IRAC Color— Color plot for asteroids and other sources 
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Fig. 4. — IRAC colors of astronomical objects and asteroids. The plus signs correspond to 
asteroids in the FLS 5° field, asterisks correspond to asteroids in the FLS 0° field and boxes 
correspond to asteroids in the SWIRE fie ld. This plot is an adaptation of lLacy et al.l (120041 ) 
as modified by iRyan fc Woodward! (120061 ). 



